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Howes et al. reply: Matthaeus, Servidio, and 
Dmitruk (hereafter MSD) [l[ comment on our Letter 
arguing that our resolution is insufficient to claim (i) that 
our gyrokinetic (GK) simulations capture simultaneously 
the MHD and kinetic sub-ion-gyroscalc plasma turbu- 
lent cascades. Our spatial resolution (64 x 64 x 128) 
is insufficient to determine converged values of spectral 
indices — a claim we do not make — but is sufficient to 
capture the nonlinear energy transfer from large to small 
scales, see a qualitative transition of electric and mag- 
netic energy spectra at the ion gyroscale, and yield re- 
sults consistent with theory above and below this scale 
(see Figs. 2 and 3 of 0). MSD speculate that the low- 
wave-number spectrum in our simulations is imposed by 
the driving. The driving occurs at the lowest mode num- 
bers (k x ,k y )L± /(2?r) = (1,0), (0,1), and (-1,0); higher- 
wave-number amplitudes are determined by nonlinear in- 
teractions. 

Our Letter presents the first kinetic simulation show- 
ing that the spectra of anisotropic turbulence in a weakly 
collisional plasma are qualitatively consistent with mea- 
surements in the solar wind (SW) 0] , a result not a pri- 
ori obvious (Q addresses a different problem). This first 
step in tackling SW turbulence at kinetic scales suggests 
that GK modeling is a promising line of inquiry. 

MSD allege that we claim (ii) that dissipation is re- 
quired to achieve our result. This apparently arises from 
a misunderstanding of the nature of kinetic Alfven waves 
(KAW), not a dissipative effect but dispersive waves 
weakly Landau-damped at k±pi > 1. We emphasize 
that damping is small for the parameters used in our 
simulation (see Fig. 1 of Q), so the spectrum attributed 
to KAW turbulence is near-dissipationless (agreeing with 
suggestions in Refs. [1,0], which MSD incorrectly claim 
to contradict our results). 

Although GK cannot describe all features of SW tur- 
bulence (as we acknowledge Q), it is a rigorous approx- 
imation for small-scale anisotropic low-frequency fluctu- 
ations (see [3] and references therein). The anisotropy 
assumption is based on a wealth of theoretical, numeri- 
cal, and observational evidence (see [2|). 

While questioning the validity of GK, MSD appear to 
suggest that Hall MHD (HMHD) is a preferable model for 
SW turbulence. Much has been learned from this tempt- 
ingly simple model [8j. The fully nonlinear anisotropic 
limit (fc|| <C k±) of HMHD can be derived from GK in 
the limit $ <§; 1 and f3 e ~ 1 (see Appendix E of Ref. [§]). 
The transition from the MHD limit to KAW occurs at 
h ~ d%/ *Jl + 2//3 e , where di — pi/y/Pl is the ion inertial 
scale. For l t 3> Pi (enabling a fluid approximation), we 
must have /9j < 1 and Tj <C T e — it is well known from 
kinetic theory fioj that this cold ion limit is the only one 
in which HMHD is strictly valid, a limit not universally 


applicable to the SW. Thus, neither the transition to the 
dispersive waves nor the effects of collisionless damping 
(not captured by any fluid-like limit) are correctly de- 
scribed by HMHD in the SW parameter regime. 

While it is true that the transition to dispersive waves 
in HMHD causes flatter electric and steeper magnetic en- 
ergy spectra, it is quite unclear whether the simulations 
MSD present are directly applicable to the SW (much 
more unclear, we believe, than in the case of GK). With- 
out a kinetic simulation to assess the importance of finite- 
ion-gyroscale effects and collisionless damping, a numer- 
ical HMHD simulation may be inconsistent with the un- 
derlying kinetic plasma physics, little more than an ad 
hoc model. It is possible to get what appears to be the 
right answer for the wrong reason, as demonstrated by 
the 1-D HMHD result shown by MSD: surely they are 
not claiming that SW turbulence is accurately modeled 
as a 1-D cascade in fc|| ! That qualitatively correct spec- 
tra seem to emerge in such studies does not invalidate 
the much more general GK approach or remove the need 
for kinetic simulations. 

Obviously, the success of GK in qualitatively repro- 
ducing the observed spectra does not in itself prove 
its applicability in space. We make no such claim and 
accept that further analytical and numerical work is 
needed. In presenting the first numerical results of our 
long-term project on kinetic turbulence, we phrased our 
conclusions carefully. The hyperbolic claims attributed 
to us by MSD — e.g., that we prove the applicability of 
GK to the SW, or show that the spectra reported in 
Ref. [H require dissipation — are an incorrect representa- 
tion of our work. 
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